Scaling of the conductance in a quantum wire 
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The conductance G of an interacting nano-wire containing an impurity and coupled to non- 
interacting semi-infinite leads is studied using a functional renormalization group method. We 
obtain results for microscopic lattice models without any further idealizations. For an interaction 
which is turned on smoothly at the contacts we show that one-parameter scaling of G holds. If 
abrupt contacts are included we find power-law suppression of G with an exponent which is twice as 
large as the one obtained for smooth contacts and no one-parameter scaling. Our results show excel- 
lent agreement with the analytically known scaling function at Luttinger liquid parameter K = 1/2 
and numerical density-matrix renormalization group data. 
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The interplay of electron correlations and a single im- 
purity in one-dimensional electron systems leads to strik- 
ing effects in the low-energy physics. 1-5 Using bosoniza- 
tion and including a single impurity in an idealized man- 
ner the application of a perturbative (in the strength of 
the impurity) renormalization group (RG) method to the 
resulting local sine-Gordon model (LSGM) led to a sim- 
ple picture: at low energy scales physical observables be- 
have as if the system is split in two chains with open 
boundary conditions at the end points. 5 The bulk part 
of the model studied is known to capture the universal 
bulk Luttinger liquid (LL) physics, 6 which is character- 
ized by the interaction dependent LL parameter K < 1 
(repulsive interaction). Considering the conductance G 
as a function of the temperature T and the strength of 
the impurity v it was argued that for fixed K the RG 
flow from weak to strong impurity strength determines a 
scaling function Gk{x) on which the data for different T 
and v can be collapsed (one-parameter scaling). Going 
beyond the perturbative RG it was shown that scaling 
indeed holds for the LSGM. 5,7 ' 8 Applying the thermody- 
namic Bethe ansatz 8 (BA) Gk was determined explic- 
itly for K = 1/2 and K = l/3. 5 < 7 ' 8 In Refs. 5, 7, 8 the 
transport in fractional quantum Hall effect systems was 
studied and thus an infinite LL without leads was con- 
sidered. Then in the impurity free case the conductance 
is renormalized by the interaction, G* = Ke 2 /h. If the 
connection of a finite LL quantum wire to non-interacting 
leads is modeled by a position dependent LL parameter 
K(x) which takes the value one in the leads (local LL 
description), the conductance is given by G* — e 2 /h in 
the impurity free case. 9,10 The question how the one- 
parameter scaling of the conductance of the system with 
an impurity is modified in the presence of leads was not 
addressed before. 

Using numerical methods it was shown that the LSGM 
captures the behavior observed in a microscopic lattice 
model in the limits of weak and strong impurities. 11 
Based on these results it is generally believed that (i) 
one-parameter scaling holds for a large class of models 



of correlated electrons with a single impurity and (ii) 
the data always scale on the universal scaling function 
Gk found within the LSGM. The issue of universality is 
of special importance since transport experiments have 
been interpreted in terms of scaling. 12 

Here we study the transport through an interacting 
nano-wire with an impurity connected to non-interacting 
semi-infinite leads. The wire is modeled by the lattice 
model of spinless fermions with nearest-neighbor inter- 
action and a hopping impurity. Using a fcrmionic func- 
tional RG method 13 we are able to calculate G. We show 
that for a single impurity, if the interaction is turned 
on very smoothly starting at the contacts and no one- 
particle scattering terms (see Eq. (1)) at the contacts 
are considered (smooth contacts) one-parameter scaling 
holds also in the presence of leads. For K =1/2 the data 
perfectly match the LSGM curve Gk=i/2 provided the 
impurity free conductance G* is taken to be e 2 /h instead 
of Ke 2 /h 5 as appropriate for systems with leads. The 
latter gives us confidence that, despite the approximate 
treatment of the interaction our approach captures the 
relevant physics. This is further supported by a com- 
parison to numerical data for the conductance obtained 
by the density-matrix renormalization group (DMRG) 
method. In the more generic case with abrupt contacts, 
i.e. if the interaction is turned on more rapidly or one- 
particle scattering terms at the contacts are included, 
we find a low-energy power-law suppression of G with 
an exponent which is twice as large as the one obtained 
for smooth contacts and no one-parameter scaling. Evi- 
dently in most experiments with quantum wires the con- 
tacts are abrupt. 

The model we study is given by the Hamiltonian 

oo 

H = - E'(sVi+ h - c -) 

j=~oo 

-M {c\c x + h.c.) - t r (ch w c Nw+1 + h.c.) 

~ P ( C AW2 C AW2+1 + h ' C -) 



1 



in standard second-quantized notation. The sum in the 
first line runs over all lattice sites j excluding j = 0, 
Nw/2, and iVw, which is indicated by the prime. The 
hopping matrix element and the lattice constant are set 
to 1. We consider an even number Nw of lattice sites 
in the interacting region (the wire). A single hopping 
impurity of strength p < 1 and modified hopping matrix 
elements ti/ r < 1 in and out of the wire are included. 
Here we focus on the half-filled system. The nearest- 
neighbor interaction Uj close to the contacts at sites 1 
and Nw is assumed to be spatially varying. The constant 
bulk value of the interaction is denoted by U. The bulk 
LL parameter is given by K = [| arccos (— -j)] , for 
U < 2, as follows from the BA solution. 14 ' 15 

In linear response and at T = the conductance can 
be determined from the one-particle Green function of 
the interacting wire taken at the chemical potential and 
calculated in the presence of the non-interacting semi- 
infinite leads. 16 For the half- filled case we obtain G = 
Ae 2 tft 2 |GWiv,i(0)| jh, with the Green function Gn w ,i- 

In our earlier applications of the RG method 13 ' 17 we 
neglected the flow of the two-particle vertex and consid- 
ered the flow of the self-energy only. Within this ap- 
proximation the LL exponents relevant for the impurity 
problem turned out to be correct to leading order in the 
interaction. We here go beyond this and include the flow 
of the vertex in an approximate way. We replace the 
three-particle vertex by its initial value 0. The flow equa- 
tion for the two-particle vertex then reads 

d A r A (a,0 :1 ,S) = 

-Tr{p A r A (...,... ;l ,S) [G A ] ' T A (a, /?;...,.. .)} 
-Tr{p A r A ( a ,...; 7 ,...)G A F A (/3,...;5,...) 

— [a <-> 0\ — [7 <-> S] + [a <-> 0, 7 «-> 5] | . 

The Greek letters stand for the quantum numbers of the 
basis in which the problem is considered and the Matsub- 
ara frequencies. On the right hand side T A is understood 
as a matrix in the variables which are not written. P A = 

G A (d A [G°' A ] -1 ) G A , with G A = ([G°' A ] _1 - £ A ) _1 

and the cut-off dependent self-energy S A . G°' A is the 
non-interacting impurity free propagator supplemented 
by an infrared cut-off. As before 13 we use a frequency 
cut-off G°' A = @{\ui\- A)G°, with A e]oo,0]. In the flow 
of the vertex we replace G A by G 0,A . We neglect the 
frequency dependence of the vertex which leads to a fre- 
quency independent self-energy. Due to this the bulk LL 
properties of the model are only partially captured by our 
approximation. In particular we miss the bulk anoma- 
lous dimension which is small compared to the impurity 
contribution included in our approach. We parameter- 
ize T A by an effective nearest-neighbor interaction with 



a renormalized amplitude U A , whose flow is determined 
by projecting onto the Fermi points. In this way the 
fixed point coupling is guaranteed to be correct to or- 
der U 2 . The coupling to the leads is neglected. Within 
these approximations the flow equation for U A closes. In 
the thermodynamic limit it can be integrated analytically 
leading to 

U A = ( AU_U2+^y\ 
U \ 2tt 27TV4TAV V ' 

Details of this approximation scheme are presented 
elsewhere. 18 The one-particle Green function of the in- 
teracting wire in the presence of the semi-infinite leads 
and thus the conductance can then be determined numer- 
ically by integrating Eqs. (14) and (15) of Ref. 13 using 
U A instead of U. On the right hand side of Eq. (14) the 
self-energy on sites 1 and Nw has to be modified: due to 
the coupling to the leads terms which can be expressed 
by ti/ r and the ti/ r = Green function of the leads at 
sites and Nw + 1 have to be added. The initial con- 
dition for S A is given by the impurity potential. Below 
we show that the above approximation scheme leads to 
exceptionally good results. 




FIG. 1. Conductance as a function of U for Nw = 12, 
p — 0.5, tu r = 1. The interaction is turned on sharply. 

It was very recently suggested that for interacting sys- 
tems the conductance through a wire can be determined 
from the persistent current which is observed in the pres- 
ence of a magnetic flux piercing a ring in which the wire 
is embedded. 19-21 In Ref. 17 two of us presented results 
which strongly support this prediction. This enables us 
to obtain the conductance using the DMRG method (for 
details see 17 and 20 ) and by comparison to verify the qual- 
ity of our approximate treatment of the interaction. In 
Fig. 1 we present RG (solid line) and DMRG 22 (circles) 
results for G/(e 2 /h) as a function of U for Uj = U, 
j = 1,2, ... , Nw — 1, h/ r = 1) an d p = 0.5. The good 
agreement of the data for interactions as large as U = 2 
gives us confidence that the functional RG with the ap- 
proximations discussed above provides data for G which 
are very close to the exact value. To determine U A we 
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numerically integrated the finite size flow equation for 
the vertex, but even for N\y — 12 the difference between 
using this and the N\y — » oo result of Eq. (2) is only 
marginal. 




FIG. 2. Conductance as a function of Nw for U = 1.5, 
p = 1, ij/ r = 1 (impurity free; circles), p — 0.5, ti/ r = 1 
(squares), and p = 0.5, ti/ r = 0.7 (diamonds). The interaction 
is turned on sharply. 

Assuming an interaction which is turned on (and off) 
sharply (in space), i.e. Uj = U for j = 1,2,..., N\y — 1, 
leads already in the absence of single particle scattering 
terms (i.e. ti/ r = p = 1) to a conductance which depends 
on U and An/. 23 ' 20 ' 17 For large Nw, i.e. at low energy 
scales, G goes to zero with a power-law as shown in Fig. 2 
(circles). This strong influence of an abrupt contact has 
not been discussed before. To unambiguously determine 
the exponent exceptionally large wires of up to 10 6 lattice 
sites have to be considered. 24 For smaller U even larger 
Nw are required. We find that G <~ N W 1_1 ^ K \ Since 
our approach involves an approximate treatment of the 
interaction we only obtain an approximation Arc for K. 
For U = 1.5 we find Arg = 0.643, in excellent agreement 
with the BA result K — 0.649. The above exponent is ex- 
actly the one found in the LSGM including one impurity 5 
if T is replaced by 1 /Nw (see below) . The other curves 
of Fig. 2 will be discussed later. To avoid the suppression 
of G due to the contacts and to investigate the role of a 
single impurity in the wire we now turn on the interac- 
tion smoothly: Uj = U{ir/2 + arctan(s[j — j s /2])]/n for 
j = 1,2,..., ] s - 1, U 3 = U for j = j s , . . . , N w /2, and 
u N w -j = Uj for j = N w /2 + 1, . . . , N w - 1. The larger 
Nw and U the smoother Uj has to be varied to obtain the 
impurity free conductance for ti/ r = p = 1. This proce- 
dure enables us to quantify the terms "perfectly" and 
"adiabatically" connected used in the field theoretical 
modeling of transport through an impurity free LL. 9 ' 10 
In the absence of an impurity we find G* = e 2 /h in agree- 
ment with the local LL description mentioned earlier. 9 ' 10 
For the results presented below we choose s and j s such 



that for ti/ r = p = 1 the relative deviation of the con- 
ductance from e 2 /h is less than 2 x 10~ 4 . For the system 
sizes considered in the following (Nw > 256) we do not 
expect the local change of the interaction strength over 
a few lattice sites to have a relevant effect on the flow 
of the bulk two-particle vertex and thus neglect this lo- 
cal change in the flow of the vertex. It turned out that 
as long as the switching on of the interaction is smooth 
enough and the bulk part of the wire is large compared 
to the switching region, G is independent of the details of 
the switching procedure as expected. To determine U A 
we use Eq. (2) obtained for Nw — * oo. 

In the LSGM the conductance is studied as a function 
of temperature. 5 ' 7 ' 8 Here we limit ourselves to T = 
but treat wires of finite length. We expect that the tem- 
perature scaling can directly be translated into a scal- 
ing in 1/Nw and thus study G as a function of l/Nw- 
For a certain value of U we later confirm this by a di- 
rect comparison of the scaling function calculated by us 
with Gk=i/2 obtained within the LSGM. The variable in 
which scaling is expected is Nw/N n , where N denotes a 
non-universal length scale, i.e. No depends on the details 
of the model and its parameters [here = N (U, p)]. 
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FIG. 3. Scaled conductance for an interaction which is 
turned on smoothly with the bulk strength U = 1 and 
ti/ r = 1. Different symbols stand for different p and 
N w = 256, 512, . . . , 4096 in each case. 

For the universality of scaling to hold it is necessary 
that G for small impurity strength, i.e. 1 — p <C 1, and 
small Nw scales as 1 — G/(e 2 /h) <~ N^ 1 ^^ and for 
large impurity strength, i.e. /) < 1, and large Nw as 
G/(e 2 /h) - N w ^ 1/K) . 5 Using the functional RG we find 
power-law behavior in both these limits with exponents 
which according to the above two relations can be ex- 
pressed consistently in terms of a single approximate LL 
parameter Krq. For example we find K^q(U — 0.5) = 
0.858 and K KG (U = 1) = 0.741 both in excellent agree- 
ment with the exact values K = 0.861 and K = 0.75 
obtained from the BA. To demonstrate one-parameter 
scaling we calculated G for p = 0.99, 0.98, 0.95, 
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p = 0.9, 0.85 . . . , 0.05, N w = 256, 512, . . . , 4096, j s = 22, 
s = 2, and various U. As an example Fig. 3 shows that 
for U = 1 the data can be collapsed on a single curve 
G RG /{e 2 /h) with x = [Nw/NoiU^)} 1 -^ as a scaling 
variable. The limiting behavior discussed above leads to 
the asymptotic x dependence (solid lines) indicated in 
the figure. 



reached. Using a site impurity, a combination of site and 
hopping impurities, and different positions of the local 
impurity we have verified that scaling holds for generic 
types of impurities and that independently of the type of 
impurity for fixed U the same scaling function is found. 
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x=[N w /N (U,p)]' 



FIG. 4. The same as in Fig. 3 but for U = 2.23. The solid 
line is the K — 1/2 scaling function of the LSGM. 



After showing that we find one-parameter scaling for 
the lattice model considered also in the presence of leads 
(provided the contacts are assumed to be smooth) we 
now compare the resulting scaling function to the one 
determined analytically within the LSGM for K = 1/2. 5 
We find that U = 2.23 leads to K RG w 1/2. This is 
fairly close to U — 2 which in the BA solution gives 
K = 1/2. Within our approximation the RG method 
does not capture the charge density wave ordering tran- 
sition at U — 2 occurring in the exact treatment of the 
impurity free model. 15 A comparison of the RG data (p, 
Nw, j s , and s as above) with the LSGM function (solid 
line) is presented in Fig. 4. The excellent agreement 
shows that in both physical situations the same scaling 
function is found if the conductance is divided by e 2 /h in- 
stead of Ke 2 /h, which is appropriate when no leads are 
present. The assumed equivalence of T and 1/N scal- 
ing holds and our method captures the relevant physics 
quantitatively despite our approximate treatment of the 
interaction. Obviously (see Figs. 3 and 4) the scaling 
function depends on the interaction, i.e. the LL param- 
eter, as expected from the LSGM. This has to be con- 
trasted to the fermionic RG procedure of Ref. 25, where 
the scaling function for arbitrary interaction turns out 
to be the non-interacting one, if the scaling variable x is 
defined as above. 

For U = 0.5, 1, and 2.23 the length scale N is shown 
in Fig. 5 as a function of the non-interacting reflection 
amplitude \R\ = (1 - p 2 )/(l + p 2 ). For A > 1 this 
scale provides a measure for how large Nw has to be 
for a given p and U before the strong impurity limit is 
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FIG. 5. No as a function of \R\ for different U. 

We next study the more generic case of a single im- 
purity in an interacting wire taking abrupt contacts into 
account. They are modeled considering either ti/ r < 1 or 
an interaction which is turned on sharply, i.e. Uj = U for 
j = 1,2, ... , Nw — I- For the latter case G as a function 
of Nw is shown in Fig. 2 (squares) for p = 0.5, tn r = 1, 
and U = 1.5. We find that the suppression of G in the 
low energy limit again follows a power-law but with the 
exponent 4(1 — 1/K) which is twice as large as the one 
obtained for smooth contacts. The same holds if the 
contacts are modeled by additional reduced t\ i r indepen- 
dent of whether the interaction is turned on smoothly or 
sharply. In this case the asymptotic regime is reached for 
much smaller systems as can be seen in Fig. 2 for p = 0.5 
and ti/ r = 0.7 (diamonds). In both the above cases the 
data for different p and N w cannot be collapsed on a 
single curve by a one-parameter scaling ansatz. Consid- 
ering other types of impurities, asymmetric coupling to 
the leads, and different positions of the local impurity 
we have verified that the exponent 4(1 — 1/K) is found 
generically. We thus believe that the low-energy and reso- 
nant tunneling 5 properties in most transport experiments 
on interacting quantum wires including a single impurity 
and coupled to non-interacting leads via contacts are gov- 
erned by the exponent 4(1 — 1/K) instead of 2(1 — 1/K). 
One-parameter scaling cannot be expected to hold. This 
observation might be relevant for the attempts to in- 
terpret recent resonant tunneling experiments on carbon 
nanotubes. 26 

In summary, using a functional RG method we 
have shown that in a lattice model of interacting 
one-dimensional electrons with a single impurity one- 
parameter scaling of the conductance holds even in 
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the presence of non-interacting leads if smooth contacts 
are considered, that is if the interaction is turned on 
smoothly and no additional one-particle scattering terms 
at the contacts are included. By comparison to numeri- 
cal DMRG data and the K = 1/2 scaling function of the 
LSGM we have verified that the approximate treatment 
of the interaction is valid up to fairly large interactions. 
We have shown that in the impurity free case a sharp 
onset of the interaction leads to a power-law suppression 
of G with an exponent 2(1 — 1/K) which is the same as 
the one found for a single impurity within the LSGM. 
In the generic case of abrupt contacts and a single im- 
purity in the wire we do not find one-parameter scaling 
and the low-energy properties are governed by the expo- 
nent 4(1 — 1/K). The functional RG approach has the 
important advantage that it is very flexible and allows to 
determine the conductance for microscopic lattice models 
of correlated electrons without any further idealizations. 
The spin degree of freedom can be included. 18 For further 
studies the combination of the fcrmionic lattice descrip- 
tion and the functional RG approach allows for a more 
realistic microscopic modeling of contacts and leads. 



ACKNOWLEDGMENTS 



15 This implies that the smallest K obtainable within the 
present model and for half-filling is K = 1/2. 

16 J.S. Langer, Phys. Rev. 127, 5 (1962); A. Oguri, Phys. 
Rev. B 56, 13422 (1997). 

17 V. Meden and U. Schollwock, Phys. Rev. B 67, 193303 
(2003). 

18 S. Andergassen, T. Enss, V. Meden, W. Metzner, U. 
Schollwock, and K. Schonhammer, to be published 

19 O.P. Sushkov, Phys. Rev. B 64, 155319 (2001). 

20 RA. Molina et al, cond-mat/0209552. 

21 T. Rejec and A. Ramsak, cond-mat/0301608. 

22 To obtain the conductance using DMRG non-interacting 
leads have to be added and an extrapolation in the length 
of the leads is required. Therefore the size of the interacting 
wire is fairly limited. 17 ' 20 

23 A. Oguri, Phys. Rev. B 59, 12240 (1999); tbid. 63, 115305 
(2001). 

24 Due to an optimization of our code suggested to us by T. 
Enss we are now able to treat much larger chains compared 
to our earlier studies 13 . E.g. obtaining G for Nw = 2 20 
takes roughly 15 minutes on a 2.4 GHz Pentium 4. We also 
thank X. Barnabe-Theriault for discussions. 

25 K.A. Matveev, D. Yue, and L.I. Glazman, Phys. Rev. Lett. 
71, 3351 (1993). 

26 H. Postma et al, Science293, 76 (2001). 



We thank A. Rosch for urging us to clarify whether 
the expected scaling behavior of the conductance is 
captured by our RG approach. U.S. is grateful to 
the Deutsche Forschungsgemeinschaft and Die Junge 
Akademie and V.M. to the Bundcsministerium fur Bil- 
dung und Forschung for support. 



1 A. Luther and I. Peschel, Phys. Rev. B 9, 2911 (1974). 

2 D.C. Mattis, J. Math. Phys. 15, 609 (1974). 

3 W. Apel and T.M. Rice, Phys. Rev. B 26, 7063 (1982). 

4 T. Giamarchi and H.J. Schulz, Phys. Rev. B 37, 325 (1988). 

5 CL. Kane and M.P.A. Fisher, Phys. Rev. Lett. 68, 1220 
(1992); Phys. Rev. B 46, 15233 (1992). 

6 F.D.M. Haldane, J. Phys. C 14, 2585 (1981). 

7 K. Moon et al, Phys. Rev. Lett. 71, 4381 (1993). 

8 P. Fendley, AWW. Ludwig, and H. Saleur, Phys. Rev. 
Lett. 74, 3005 (1995). 

9 I. San and H.J. Schulz, Phys. Rev. B 52, R17040 (1995). 

10 D.L. Maslov and M. Stone, Phys. Rev. B 52, R5539 (1995); 
D.L. Maslov, Phys. Rev. B 52, R14368 (1995). 

11 S. Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992). 

12 F.P. Milliken, OP. Umbach, and RA. Webb, Solid State 
Commun. 97, 309 (1996); M. Bockrath et al, Nature 397, 
598 (1999); Z. Yao et al, Nature 402, 273 (1999). 

13 V. Meden, W. Metzner, U. Schollwock, and K. Schon- 
hammer, J. of Low Temp. Physics 126, 1147 (2002). 

14 F.D.M. Haldane, Phys. Rev. Lett. 45, 1358 (1980). 



5 



